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Description 



The invention relates to a method of analysis and 
control, in real time, of a fluid influx into a hydrocarbon 
well which occurs during drilling. When, during the 5 
drilling of a well, after passing through an imperme- 
able layer, a permeable formation is reached contain- 
ing a liquid or gaseous fluid under pressure, this fluid 
tends to flow into the well if the column of drilling fluid, 
known as drilling mud, contained in the well is notable 10 
to balance the pressure of the fluid in the aforemen- 
tioned formation. The fluid then pushes the mud up- 
wards. There is said to be a fluid influx or "kick". Such 
a phenomenon is unstable: as the fluid from the for- 
mation replaces the mud in the well, the mean density 15 
of the counter-pressure column inside the well de- 
creases and the unbalance become greater. If no 
steps are taken, the phenomenon runs away, leading 
to a blow-out. 

This influx of fluid is in most cases detected early 20 
enough to prevent the blow-out occurring, and the 
first emergency step taken is to close the well at the 
surface by means of a blow-out preventer. 

Once this valve is closed, the well is under con- 
trol, but only as long as the well pressure does not ex- 25 
ceed the formation fracture pressure, otherwise 
there can result ah underground blow-out. A choke 
valve is used at the surface to relieve, in a controlled 
manner, the pressure which has been building up in 
the well. There is a conflict between the need to close 30 
the outlet choke valve sufficiently to ensure that the 
bottomhole pressure remains high enough to be 
above the formation pressure and so avoid a further 
influx but low enough to avoid the risk of fracturing the 
formation higher up the wellbore, the result of which 35 
would be an underground blow-out. In addition the 
well pressure must build up sufficiently to be able to 
determine enough information about the influx to en- 
sure that subsequent control of the choke valve will 
be correct The information that is of particular value 40 
to the driller is: 

- The formation pressure, so that the correct 
mud weight to be used for the mud circulated 
to replace the original fluid can be selected and 

so that the choke valve can be operated to 45 
maintain downhole pressure above the forma- 
tion pressure and so ensure no further influx 
occurs. 

- Details of the influx: the crucial information is 
whether it consists of gas or water or oil. This so 
decides subsequent action in circulating out 

the influx. The density of the influx if it is gas 
and the rate at which it is rising up the annulus 
is sufficient to determine the maximum attain- 
able pressure at the casing show and so decide 55 
whether or not fracture will occur. Also the vol- 
ume of the influx is important in determining 
the subsequent well kill operations as the orig- 



inal volume estimate, which is taken from the 
surface pit gain, is notoriously inaccurate. 
- As the pressure builds up in a shut-in well, the 
influx flow rate falls away until eventually it 
ceases. It is vitally important to know when the 
influx has ceased because any further delay in 
operating the choke valve to reduce wellbore 
pressure can result in fracturing the formation. 
However a premature operation of the choke 
valve would result in a further influx of gas with 
possible disasterous consequences. 
Once the well is under control under the opera- 
tion of the choke valve, the formation fluid can be 
safely circulated out and the mud then weighted to 
enable drilling to continue without danger. If the for- 
mation fluid that has entered the well is a liquid (brine 
or hydrocarbons, for example), the circulation of this 
fluid does not present any specific problems, since 
this fluid scarcely increases in volume during its rise 
to the surface and, therefore, the hydrostatic pres- 
sure exercised by the drilling mud at the bottom of the 
well remains more or less constant. If on the other 
hand the formation fluid is gaseous, it expands on ris- 
ing and this creates a problem in that the hydrostatic 
pressure gradually decreases. To avoid fresh influxes 
of formation fluid being induced during "circulation" of 
the influx, in other words while the gas is rising to the 
surface, a pressure greater than the pressure of the 
formation has to be maintained at the bottom of the 
well. To do this, the annulus of the well, this being the 
space between the drill string and the well wall, must 
be kept at a pressure such that the bottom pressure 
is at the desired value. It is therefore very important 
for the driller to know as early as possible, during cir- 
culation of the influx, if a dangerous incident is on the 
point of occurring, such as a fresh influx of fluid or the 
commencement of mud loss due to the fracture of the 
formation. 

The usual means of analysis and control avail- 
able to the driller comprise the mud level in the mud 
tank, the mud injection pressure into the drill pipes, 
and the well annulus surface pressure. These three 
data allow the driller to calculate the volume and na- 
ture of the influx, and also the formation pressure. It 
is on this information that he bases his influx circula- 
tion programme. 

Interpreting the data nevertheless poses some 
problems. Firstly, the assessment of the volume of 
the influx, which is important in order to determine 
the nature of that influx, is inaccurate. It is in fact 
made by comparing the mud level in the tank with a 
"normal" level, ie the level that would occur in the ab- 
sence of the influx. But this reference is drff icultto de- 
termine: on one hand the mud level changes con- 
stantly during drilling, because part of the mud is 
ejected with the well cuttings; on the other, the mud 
level in the pits rises when the well is closed, because 
the mud return lines empty. The estimate of the influx 
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volume is therefore approximate. As a result, deter- 
mining the nature of the influx is also uncertain. The 
influx density calculations thus often lead to the con- 
clusion that the influx is a mixture of gas and liquid (oil 
or water) whereas it may in fact be a gas or a liquid 
only. It should also be noted that this calculation can 
not be made when the influx is in a horizontal part of 
the well. 

For all these reasons, influx analysis is not re- 
garded as a reliable technique today. 

Several methods have already been proposed for 
analysing and/or controlling fluid influxes into an oil 
well from an underground formation being drilled. For 
example, in US patent 4,867,254 the value of the 
mass of gas in the annulus is monitored in order to de- 
termine either a fresh gas entry into the annulus or a 
drilling mud loss into the formation being drilled. In EP 
patent application 0,302,558, the variations of the 
flow rate or the pressure of the inlet drilling mud are 
compared with the variations of the flow rate or the 
pressure of the outlet mud and, from the comparison, 
the nature and volume of the influx are determined. 
Other examples of methods for detecting and/or con- 
trolling a fluid influx can be found in US patents 
4,840,061; 3,740,739; 3,760,891; 4,253,530 and 
4,606,415. 

However, the methods of the prior art are often 
not sufficiently accurate to allow a correct determina- 
tion of the parameters characterizing the influx and 
the well conditions. For example, the precise time to 
open the choke valve in order to control the well is 
either not described or predicted as being later than 
necessary. 

The present invention offers a method of deriving 
the required information to analyse and control a fluid 
influx in a borehole from an analysis of the surface in- 
let or outlet pressure monitored on a continuous basis 
when the well is shut-in and operating the choke 
valve at the right time and in the correct manner. The 
proposed method may be applied in deviated and 
even horizontal wells. 

More precisely, the invention relates to a method 
of real time analysis and control of a fluid influx from 
an underground formation into a wellbore being dril- 
led with a drillstring, a drilling mud circulating from the 
surface down to the bottomhole into the drillstring and 
flowing back to the surface in the annulus defined be- 
tween the wall of the wellbore and the drillstring, 
wherein the well is shut-in when the influx is detected 
and wherein the mud pressure, which is the inlet pres- 
sure Pi of the drilling mud, is repeatedly by measured 
as a function of time at surface, the method being 
characterized by determining, from the increase of 
the mud pressure measurement, the time t c corre- 
sponding to the minimum gradient in the increase of 
the mud pressure when the relationship: 

p, = A{1 - e- <V) for t ^ tc (1) 
is no longer valid and the relationship 




Pi = B + C^f-fc) for t i= tc (2) 
starts to apply, wherein A, 8, C, and C 2 are constants 
t and is the time at which time the inlet pressure p f is 
5 substantially equal to the difference between the for- 
mation pressure p f and the hydrostatic pressure p H 
the drilling mud and controlling the well from said time 

When inlet pressure p f is measured, time t c cor- 
10 responds to the time when the inlet pressure p y is sub- 
stantially equal to the difference between the forma- 
tion pressure p f and the hydrostatic pressure p H cre- 
ated by the density of the drilling mud. The formation 
pressure p f is derived by adding the inlet pressure p f 
15 at time t G to the hydrostatic pressure p H . The rate of 
change dp/ of the inlet pressure p f is monitored at time 
t^ said rate of change being compared with a prede- 
termined value, and the type of influx is determined 
from said comparison. In addition, the volume and the 
20 density of the influx can be derived from the determi- 
nation of time t^. 

The invention applies as well to analysis based 
on continuously monitored outlet pressure p D . For il- 
lustrative purposes only the inlet pressure will be 
25 mentioned from now on. 

The characteristics and advantages of the inven- 
■ tion will be seen more clearly from the description that 
follows, with reference to the attached drawings, of a 
non-limitative example of the method mentioned 
30 above. 

Figure 1 shows schematically the drilling mud cir- • 
cuit of a well during control of an influx. 

Figure 2 shows in diagram form the hydraulic cir- 
cuit of a well during control of a gas influx. 
35 Figure 3 shows an example of inlet pressure pi as 

a function of time, as predicted by the prior art and as 
observed during a numerical simulation of a gas kick, 
in accordance with the present invention. 

Figure 1 shows the mud circuit of a well 1 during 
40 a formation fluid influx control operation. The bit 2 is 
attached to the end of a drill string 3. The mud circuit 
comprises a tank 4 containing drilling mud 5, a pump 
(or several pumps) 6 sucking mud from the tank 4 
through a pipe 7 and discharging it into the well 1, 
45 through a rigid pipe 8 and flexible hose 9 connected 
to the tubular drill string 3 via a swivel 17. The mud 
escapes from the drill string when it reaches the bit 2 
and returns up the well through the annulus 10 be- 
tween the drill string and the well wall. In normal op- 
so eration the drilling mud flows through a blow-out pre- 
venter 12, which is open, into the mud tank 4 through 
a line 24 and through a vibratory screen not shown in 
the diagram to separate the cuttings from the mud. 
When a fluid influx is detected, the blow-out preven- 
55 ter 12 is closed. Having returned to the surface, the 
mud flows through a choke valve 13 and a degasser 
14 which separates the gas from the liquid. The drill- 
ing mud then returns to the tank 4 through line 15. 
The mud inlet flow rate Q| may be measured by 
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means of a flow meter 16 and the mud density is 
measured by means of a sensor 21, both of these fit- 
ted in line 8. The inlet pressure p, is measured by 
means of a sensor 18 on rigid line 8. The outlet pres- 5 
sure p 0 is measured by means of a sensor 1 9 fitted be- 
tween the blow-out preventer 12 and the choke 13. 
The mud level in the tank 4 is measured by means of 
a level sensor 20 fitted in the tank 4. This level will in- 
crease if a kick is taken and this pit gain is a simple 10 
and basic estimate of the volume of the influx. The 
sensors are connected to a data acquisition and proc- 
essing system 22. 

In order to exploit the present invention it is suf- 
ficient to measure at least pi or p 0 during the shut-in is 
phase, before the choke valve begins to be operated. 

Figure 2 represents in simplified form hydraulic 
circuit of a well when the operator is preparing to cir- 
culate the fluid influx 30 that has entered the well. 
The gas influx 30 produced by the formation being 20 
drilled has been represented, rising in theannulus 10. 
The arrows represented in the drill string 3, the drill 
bit 2 and the annulus 1 0 indicate the circulation of the 
mud when the pumps 6 are working. Immediately af- 
ter detecting an influx, the pumps are shut down and 25 
the blow-out preventer 12 and choke 13 are closed. 
The well is thus isolated or "shut-in" and the drilling 
mud is immobolized in the well. The driller then meas- 
ures at the surface the inlet pressure pj in the pipes 
by means of the sensor 8 and the outlet pressure p Q 30 
in the annulus by means of sensor 19 between the 
wellhead and the control choke 13. 

For the sake of clarity in explaining the method 
it is assumed here that the section of the annulus has 
a constant area A from the bottom to the top of the 35 
well. But the method may be used even if this section 
is not of constant area. 

Figure 3 shows the variations of the mud inlet 
pressure p, as a function of time t during a kick, which 
is detected and controlled by closing the blow-out pre- 40 
venter 12. Curve 30 represented in plain line repre- 
sents the usual field expectation of variation of p h and 
curve 32 in dashed lines represents the expected va- 
riation of pi in accordance with the present invention. 
The kick begins at time t v Before that the inlet pres- 45 
sure pi is relatively constant From time to t 2 , Pi de- 
creases very slightly until time t 2 when the kick is de- 
tected. The period of time (t 2 - U) between the start 
of the kick and its detection could be, say 5 minutes 
depending on formation productivity. At time t 2 , the so 
mud pumps are stopped. The inlet pressure p, falls 
sharply during a few seconds down to a minimum 
pressure p m m at time t 3 . The blow-out preventer is fully 
closed at time t» which is usually called the shut-in 
time. The elapsed time between t 2 and t> is the time 55 
it takes to close the blow-out preventer (about 1 min- 
ute). At time t, the inlet pressure p, rises until it reach- 
es a constant value equal to the difference between 
the formation pressure Of and the hydrostatic pres- 



sure p H . The period oTEime to reach this value is of the 
order of 5 to 10 minutes depending on formation pro- 
ductivity and includes the recovery time of the forma- 
tion. The well is shut-in completely since the pumps 
6 are stopped and the blow-out preventer 12 and the 
choke valve 13 are closed. From the time tg the well 
is shut-in, the pressure p, begins to increase for two 
reasons: 

a) The mass of the fluid influx in the wellbore 
keeps increasing as long as more and more fluid 
is produced by the formation into the wellbore. 
Since the volume of the wellbore is constant, the 
pressure p ( will increase until the influx shuts it- 
self off. 

b) If the influx is gas, it rises up the annulus at 
some slip velocity relative to the mud. As it rises 
within a fixed volume (the well is shut-in), the 
pressure increases as the gas can only expand 
a very limited amount. 

The manner in which the pressure builds up is a 
function of the volume and compressibility of the mud 
system and of the influx, the rate at which the influx 
was flowing from the formation when the blow-out 
preventer was closed as well as the rate of rise of the 
influx fluid in the annulus if it is gas. 

It is usual field practice, in recognition of phenom- 
enon (a) above, to wait until the surface pressure 
ceases to increase (when pi = p r p H on Figure 3 after 
the time y and to identify this instant as the time at 
which the influx ceased. From the value of the sur- 
face pressure at this time, the formation pressure, the 
influx density and the manner in which to control the 
choke valve are determined. 

However all of this information is deficient in the 
case of a gas influx, as recognised in the present in- 
vention, because the shut-in pressure never actually 
ceases to increase due to the phenomenon (b) above 
mentioned. The influx density calculation can conse- 
quently be grossly in error and the formation pressure 
estimate wrong. 

The usual field practice described above stems 
from the knowledge that the bottomhole pressure p* 
is lower than the formation pressure Df (since an in- 
flux is flowing from the formation into the borehole) 
and the bottomhole pressure Pw increases until it 
meets the formation pressure p f , beyond which time 
there is no further influx of fluid from the formation 
into the borehole. At that time, the inlet pressure p, is 
equal to the formation pressure Df minus the hydro- 
static pressure p H . The formation pressure p, and the 
hydrostatic pressure p H being constant, the inlet pres- 
sure pj reaches a constant value equal to (Df - p H . This 
is illustrated by the curve 30 in plain line, after the 
time tg, on Figure 3. 

However, this is not realistic and the inventor of 
the present invention has demonstrated that in fact 
the inlet pressure p, can be given by the following two 
equations: 
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p, = A(1 - e-*) (1) 
from the time t = ^ to the time t = tc, and 
P. - B + c, (t - y (2) 
from and after the time tc. wherein A, B, q and C2 are 
constants. 
By taking: 

A = (Pf - p H ) + Ci/Ca 

B = (Pr - p H ) 
the two first equations become: 

Pi = {(Pf - Ph) + ~K1 • e-*} (1a) 

from the time t = tg to the time t = tc, and 

Pi - (Pt - Ph) + c, (t - y (2a) 
from and after the time tc. 

The time tc is defined as the time when the influx 
stops and therefore the time when Pf = Pw and C2 
are constants defined hereafter. 

A non-uniform geometry modifies the detail of 
these expressions but not the principle being descri- 
bed. 

As a fact, it is then necessary to add on the right 
member of equation (2) a third team equal to + E (t - 
y 2 , E being an arbitrary coefficient introduced to ac- 
count for the departure from linearity of equation (2) 
caused by changes in area as the gas leaves the re- 
gion of the drill collars. 

In Figure 3, curve 32 in dash lines represents the 
variation of inlet pressure pi during a shut-in period, in 
accordance with equations (1) and (2). 

The time ^ can be determined directly from the 
measurement of the inlet pressure p jt as the inflection 
point 34 of curve 32 or the point of minimum gradient. 
This is because the minimum gradient in the increase 
of ^ versus time occurs precisely for t = tc (point 34). 
The determination of the minimum gradient can be 
done for example by plotting the curve 32 with the 
pressure measurement versus time or with the help 
of a computer. 

Another way to determine ^ is to do it by compu- 
tational means. The way to do so is to match the 
measured data pi versus time with predictions of p, 
from equations (1) and (2) based on assumed values 
of c lP C2, A and B and refining the assumed values un- 
til a good match is obtained. The match is obtained 
when the limit ^ is found for the two equations (1) and 
(2), when equation (1) is not valid anymore and equa- 
tion (2) starts to apply. The same curve fitting process 
can obviously apply starting from equations (1a) and 
(2a). 

When equation (1) applies, for times less than tc, 
we have unknown parameters p f , c 1( C2 and when 
equation (2) applies, for times exceeding tc, we have 
unknown parameters c 1 and Pf. 

A value for the time ^ is first assumed. Then it is 
a straightforward matter to use least squares or some 
other appropriate curve fitting method to determine 



Pf, c 1t C2 from the region o<t<tc, comparing measure- 
ments with predictions of equations (1a) and c 1f Pf 
from the region tc<t comparing measurements with 
5 predictions of equation (2a). Having done this with 
the assumed value of tc, there are several further 
conditions to be met Namely that the two curves 
must coincide at time tc and that gradients of the 
curves at time ^ must match. Furthermore the para- 

10 meters Pf and c-, found from the curve fitting process 
with equation (1a), on the one hand, and from the 
curve fitting process with equation (2a) on the other 
hand, must be consistent If these conditions are not 
met, then the time ^ is adjusted. The process is re- 
ts peated iteratively until these conditions are met and 
then all the parameters ^ c 1f C2 and Pf are known. 

The determination of tc calls for several remarks. 
In the method usually applied in the field, the time t* 
to open the choke valve when p, does not increase 

20 anymore (when pj = Pf - p H on Figure 3) is difficult to 
determine with accuracy since p, rises asymptotically 
toward a plateau. The driller is therefore not really 
sure of the right instant to open the choke valve. By 
comparison, the curve 32 cuts the horizontal line (p, 

25 - p H ) at point 34, going over that line. Time tc which 
corresponds to the point of minimum gradient at the 
intersection between this horizontal line and curve 
32, is therefore easy to determine. The driller, who in 
accordance with the present invention, opens the 

30 choke valve at time tc, knows perfectly well the right 
instant to do it 

A second remark is that the inlet pressure pi con- 
tinues to increase (curve 32) after the time tc, contra- 
ry to the usual belief that Pi reaches a constant value 

35 and stops rising for a while. 

Another remark is that time tc occurs before time 
U of the prior art. As a consequence, there is a higher 
risk to fracture the formation with the usual field prac- 
tice, particularly since p t continues to increase after 

40 time tc. It is therefore very important to determine pre- 
cisely the time tc- in addition, the accuracy of other 
parameter values depends on the precision in deter- 
mining tc, since tc is used later on to determine other 
parameters. It may be noted that the use of both 

45 equations (1 ) and (2) to determine tc implies that time 
tc is passed before it is determined. This is true but is 
not of any consequence. The driller will then known 
the true ^ and would always delay for a short period 
the opening of the choke in order to give a margin of 

50 safety in controlling the downhole pressure to be not 
just at the formation pressure but marginally above it. 

When time tc has been determined, in accor- 
dance with the present invention, the inlet pressure 
Pi is determined at time tc, for example directly from 

55 the pressure measurement If no measurement was 
made at that particular time tc, then the value of p, at 
time tc is extrapolated from the measurement made 
right before and after tc- 

Then the formation pressure p f at time ^ is calcu- 
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lated in order to better control the opening of the 
choke valve and to determine the mud density suffi- 
cient to kill the well. The formation pressure is given 
by: 5 

Pf = Pi + Ph 
Any error on the determination of tc and subse- 
quently pi leads to a same error on the value of pf. This 
is important since the driller has to keep the bottom- 
hole pressure at least equal to the formation pres- 10 
sure, and therefore the inlet pressure pi large enough, 
by adjusting the opening of the choke valve. Any error 
on the value of pr leads therefore to a wrong control 
of the choke valve. The hydrostatic pressure p H is de- 
termined, as known in the art, from the density d m of 15 
the mud presently in the well and from the true vert- 
ical depth. 

It must be realized that, if the curve fitting method 
has been used to obtain tc, as explained previously, 
then the value of p f is also obtained at the same time, 
together with the values for c, and c 2 . 

In order to determine the type of influx, gas or liq- 
uid, the rate of change of inlet pressure dp, is comput- 
ed from the measured data at the time If the rate 
dpi is very small, less than say 0.03 bar/min, then the 25 
influx is not gas. This can be ascertained even in a 
horizontal well. 

In accordance with one characteristic of the in- 
vention, the volume of influx V 0 is determined. The in- 
ventor has determined that the constant q of equa- 30 
tion (1) is given by: 

= Vod m gv fl 
V G + PhX^ 
wherein p H is the hydrostatic pressure, X m the com- 
pressibility of the mud in the well, V m the volume of 35 
the mud in the well (drill string and annulus) d m is the 
density of the mud, g is the gravitational acceleration 
and v fl is the rate of rise of the gas in the annulus. The 
value of v Q is obtained from experimental conditions 
in flow simulators and is therefore known. *o 

From the time derivative of equation (2), the rate 
of change dpj of inlet pressure is: 
dp, = d 

Ct can thus be determined by determining the 
rate of change of pj at time 45 

Writing equation (3) for V G and substituting c-j by 
dpi, one obtains: 

v PHX m V m dp, 
d m gv fl -dp, 

The compressibility X m of the mud is known or 50 
can be determined easily. The rate of rise dpj of the 
inlet pressure has been determined previously and, 
the other parameters of equation (4) being known, 
then the value of V 0 can be computed. The volume of 
influx so determined is a better estimate than the one 55 
obtained with the usual pit gain measurement. 

There will be situations where the operator is 
more confident in the pit gain measurement than in 
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2rivea from 



the value of V 0 deriveo from equation (4) wherein an 
estimate of rate of rise of gas v g is inferred. In that 
case, the value of v B is obtained from equation (4) us- 
ing for V 0 the pit gain. 

However, if the difference (Q D - Q,) between the 
outlet flow rate Qo and inlet flow rate Q { has been 
measured between times t 2 and tg on Figure 3, then 
the volume V 0 of influx can be estimated from the fol- 
lowing expression, derived by the inventor, of the con- 
stant c 2 of equation (1): 

Qi) 1 



(p,- Ph ) Vo/Ph + 



(5) 



so that: 



V 0 = p H 



(Qq - Qi) 



PhX^ (6) 



C2 (Pf - Ph) 

According to a further aspect of the invention, the 
density of the influx d fl is determined, even if the well 
is deviated from the vertical, as follows: 
20 In a well with constant annulus area S the density 

of the influx d g is determined from a comparison of 
the inlet and outlet pressure at time tc. 



Po - Pi - Pfr = (d m 



d g )g-^cosa (7) 



where a is the angle of inclination of the drill collars 
from vertical, and the frictional pressure drop P fr is 
due to the relative motion of the gas with respect to 
the mud. This term is small and would be ignored if 
no test work was available to give an estimate of the 
value. 

This expression for the influx density d g will indi- 
cate whether there is gas, oil or water entering the 
well bore. When non-constant area annuli are consid- 
ered then due account would need to be taken of the 
area changes in the relationship (7). 

The prefered mode of the invention has been de- 
scribed with respect to measuring inlet pressure p,. 
However, the invention may also be practised in an 
equivalent manner by measuring outlet pressure p 0 as 
it varies with time. Equation (7) is a relationship be- 
tween outlet pressure p 0 and inlet pressure Pi during 
the shut-in, but only as long as the annulus area is 
constant This expression contains unknown terms 
such as the fluid density d 0 , the frictional pressure p fr 
and the volume of influx V Q is often poorly estimated 
by pit gain measurements. However, even with a non- 
constant annulus cross sectional area, to a good ap- 
proximation 

p 0 - pi = constant, 
where the constant is at present unknown. Thus all of 
the subsequent discussion relating to the use of p, to 
determine tc, Pf, c 1f C2 can be applied to p 0 where the 
unknown constant will be determined from the differ- 
ence 

(p 0 - p,) at time tc. 
As an example, the variation of the outlet pres- 
sure Po versus time, after the shut-in time follows 
a curve similar to curve 32 on Figure 3. From this p 0 
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curve, time tc is determined corresponding to the 
point of minimum gradient, and the values of the con- 
stants c-, and C2 are derived from this p Q curve, as be- 
fore. Then, if p, has also been measured, the forma- 
tion pressure p f , the volume V 0 and density d m are de- 
termined as previously. 



Claims w 

1. A method of real time analysis and control of a flu- 
id influx from an underground formation into a 
wellbore being drilled with a drillstring, a drilling 
mud circulating from the surface down to the bot- 
tomhole into the drillstring and flowing back to the 
surface in the annulus defined between the wall 
of the wellbore and the drillstring, wherein the 
well is shut-in when the influx is detected and 
wherein the mud pressure, which is the inlet pres- 
sure pi of the drilling mud, is repeatedly by meas- 
ured as a function of time at surface, the method 
being characterized by determining, from the in- 
crease of the mud pressure measurement, the 
time t c corresponding to the minimum gradient in 25 
the increase of the mud pressure when the rela- 
tionship: 

P/ = 4(1 . e - V) for t ^ tc (1) 
is no longer valid and the relationship 

Pi = B + d(f - y for t ^ tc (2) 
starts to apply, wherein A, 6, C, and C-i are con- 
stants and t is the time at which time the inlet 
pressure p f is substantially equal tothe difference 
between the formulation pressure p f and the hy- 
drostatic pressure p H of the drilling mud and con- 
trolling the well from said time k 

2. The method of claim 1, wherein the hydrostatic 
pressure p H is computed from the mud density d m 

and the drilled depth and wherein the formation 40 
pressure p f is derived by adding the inlet pressure 
Pi at time t c with the hydrostatic pressure p H . 

3. The method of preceding claims 1 or 2, wherein 

the rate of change d p of the mud pressure is moni- 45 
tored at time t c , said rate of change is compared 
with a predetermined value and the type of influx 
is determined from said comparison. 

4. The method of any preceding claim, wherein the so 
time tc is determined by matching the measure- 
ment of inlet pressure p f versus time with predic- 
tions of pi values from equations (1) and (2) based 

on assumed values of A, 6, C t and C 2 and refin- 
ing the assumed values until a good match is ob- 55 
tained. 

5. The method of claim 4, wherein 



35 



A = (p, - p„) + 



Ci 



B - Pf - Ph 
in which p f is the formation pressure and p H is the 
hydrostatic pressure of the mud. 

The method of any preceding claim, wherein the 
rate of change c/p/of the inlet pressure p f is deter- 
mined at time t c and the value of d is taken equal 
to said rate of change dp h 



7. The method of any one of the preceding claims, 
wherein the rate of change d p of the mud pres- 
*5 sure is determined at time t c and the volume V Q 
of the influx is computed from the equation: 
v = PHXmVmdp 

° " dngVg - dp 

in which p H is the hydrostatic pressure of the 
20 mud, X m is the compressibility of the mud, v m is 

the volume of the mud in the wellbore, d m is the 
density of the mud, g is the gravitational acceler- 
ation and v g is the mean rate of rise in the influx 
in the wellbore. 



8. The method of any one of the claims 1 to 6, 
wherein the rate of change d p of the mud pres- 
sure is determined at time t^ the pit gain volume 
V 0 is measured and the mean rate of rise v g of the 
30 influx in the wellbore is computed from the equa- 

tion: 

w _ Ph(Q 0 - Q/) Y w 
d(Pf - Ph) 

in which p H is the hydrostatic pressure of the 
mud, p f is the formation pressure, X m is the com- 
pressibility of the mud, V m is the volume of the 
mud in the wellbore. 



Patentanspruche 

1. Ein Verfahren der in Realzeit erfolgenden Analy- 
se und Steuerung eines Fluideinstroms aus einer 
untertagigen Formation in ein Bohrloch wahrend 
dessen Abteufens mittels eines Bohrstrangs, wo- 
bei eine Bohrspulung von uberTage nach unten 
zur Sohle in dem Bohrstrang zirkuliert und zur 
Oberf lache in den Ringraum, begrenzt zwischen 
der Wandung des Bohrlochs und dem Bohr- 
strang, zuruckflieBt, wobei das Bohrloch ge- 
sperrt wird, wenn der Einstrom erkannt wird, und 
wobei der Spulungsdruck, bei dem es sich urn 
den EinlaBdruck p { der Bohrspulung handelt, wie- 
derholt als eine Funktion der Zeit uber Tage ge- 
messen wird, welches Verfahren gekennzeichnet 
ist durch die Bestimmung, aus der Zunahme der 
Spulungsdruckmessungen, der Zeit % entspre- 
chend dem Minimumgradienten in der Zunahme 
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des Spulungsdruckes, wenn die Beziehung 
p, = 4(1 - e- <V) furt^tc (1) 
nicht mehr gultig ist und die Beziehung 

Pi = B + C,(t - t c ) furt^tc (2) 
wirksam zu werden beginnt, wobei A, B, G, und 
C 2 Konstanten sind und t die Zeit ist, bei welcher 
der EinlaBdruck pi im wesentlichen gleich ist der 
Differenz zwischen dem Formationsdruck p f und 
dem hydrostatischen Druck p H der Bohrspulung, 
und Steuern des Bohrlochs von der Zeit t^ 

Das Verfahren nach Anspruch 1 , bei dem der hy- 
drostatische Druck p H berechnet wird aus der 
Spulungsdichte d m und der abgeteuften Tiefe, 
und bei dem der Formationsdruck Pf abgeleitet 
wird durch Addieren des Einlafcdrucks pj zum 
Zeitpunkt tc mit dem hydrostatischen Druck p H . 
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Dichte der Spulurfij ist, g die Erdbeschleunigung 
ist und v g die mittiere Anstiegsrate des Einstroms 
in das Bohrloch ist. 

Das Verfahren nach einem der Anspruche 1 bis 
6, bei dem die Anderungsrate d p des Spulungs- 
drucks zum Zeitpunkt ^ bestimmt wird, das Tank- 
gewinnvolumen V 0 gemessen wird und die mittie- 
re Rate des Anstiegs v fl des Einstroms in das 
Bohrloch berechnet wird aus der Gleichung: 
w _ P^Qo - Q/) n v w 

worin p H der hydrostatische Druck der Spulung 
ist, Pf der Formationsdruck ist, X m die Kompres- 
sibiiitat der Spulung ist, V m das Volumen der Spu- 
lung im Bohrloch ist. 



3. Das Verfahren der vorangehenden Anspruche 1 20 
Oder 2, bei dem die Anderungsrate d p des Spu- 
lungsdrucks zum Zeitpunkt tc uberwacht wird, 
welche Anderungsrate verglichen wird mit einem 
vorbestimmten Wert und der Typ des Einstroms 
bestimmt wird aus diesem Vergleich. 25 

4. Das Verfahren nach einem der vorangehenden 
Anspruche, bei dem der Zeitpunkt ^ bestimmt 
wird durch Anpassen der Messung des EinlaR- 
drucks pj uber der Zeit an Voraussagen von p r 30 
Werten aus Gleichungen (1) und (2), basierend 

auf angenommenen Werten von A, B, C, und C2, 
und Verbessern der angenommenen Werte, bis 
eine gute Anpassung erzielt ist 

35 

5. Das Verfahren nach Anspruch 4, bei dem 

A = (Pr - p„) + £ 
B = Pf - Ph 

wobei Pf der Formationsdruck ist und p H der hy- *o 
drostatische Druck der Spulung. 

6. Das Verfahren nach einem der vorangehenden 
Anspruche, bei dem die Anderungsrate dp t des 
Einla&drucks p, zum Zeitpunkt tc bestimmt wird 45 
und der Wert von C, gleich der Anderungsrate dp, 
gewahlt wird. 

7. Das Verfahren nach einem der vorangehenden 
Anspruche, bei dem die Anderungsrate d p des so 
Spulungsdrucks zum Zeitpunkt ^ bestimmt wird 

und das Volumen V 0 des Einstroms berechnet 
wird aus der Gleichung: 

v - PHXmVmdp 

° " d^vg - dp 55 
worin p H der hydrostatische Druck der Spulung 
ist, X m die Kompressibilitat der Spulung ist, V m 
das Volumen der Spulung im Bohrloch ist, d m die 



Revendications 

1. Un proc6d6 d'analyse et de r6glage, en temps 
r§el, d'un influx de f luide venant d'une formation 
souterraine vers un puits de forage en cours de 
forage k I'aide d'une chatne de forage, une boue 
de forage circulant depuis la surface en descen- 
dant vers le fond du trou k I'inferieur de la chaine 
de forage et en s'6coulant en retour vers la sur- 
face dans I'anneau d6fini entre la paroi du trou de 
forage et la chatne de forage, dans lequel le puits 
est ferm6 lorsque I'inf lux est d§tecte et dans le- 
quel la pression de boue, qui est la pression d'ad- 
mission pi de la boue de forage, est mesurSe de 
fa$on r6p6t6e en fonction du temps, le proc6d6 
6tant caract§ris6 par les Stapes consistant k: de- 
terminer, k partirde I'augmentation de la mesure 
de pression de la boue, le temps ^ correspondant 
au gradient minimal de I'augmentation de la pres- 
sion de boue lorsque la relation: 

Pi = A(1 -e- ct )pourt=itc (1) 
n'est plus valide et que la relation 
Pi = B + C^t-y pourt^tc (2) 
commence k s'appliquer, ou A, B, d et C2 
sont des constantes et t est le temps auquel la 
pression d'admission p, est sensiblement 6gal k 
la difference entre la pression de formation Pf et 
la pression hydrostatique p H de la boue de forage, 
et k r6gler le puits k partir dudit temps t^ 

2. Le proc6d6 selon la revendication 1, dans lequel 
la pression hydrostatique p H est calcu!6e k partir 
de la densit6 d m de la boue et de la profondeur fo- 
r6e et dans lequel la pression de formation Pf est 
d6riv6een ajoutanta la pression hydrostatique p H 
la pression d'admission pi k Pinstant V 

3. Le proc6d6 selon Tune des revendications 1 ou 2, 
dans lequel la vitesse de variation d p de la pres- 
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sion de la boue est surveilleea ('instant tc, ladite 
vitesse de variation est compare k une valeur 
pr6d6termin6e et le type d'inf lux est d6termin6 k 
partir de ladite comparaison. 



4. Le proc6d6 selon Tune des revendications pr6c6- 
dentes quelconque, dans iequel le temps tc est 
d6termin§ en faisant concorder la mesure de la 
pression d'admission p { en fonction du temps 10 
avec des predictions de valeurs pi venant des 
Equations (1) et (2) bashes sur des valeurs sup- 
poses de A, B, C, et C 2 et en af f inant les valeurs 
supposes jusqu'S ce qu'une bonne concordan- 
ce so it obtenue . 15 



5. Le proc6d6 selon la revendication 4 dans Iequel 
A = (Pf-p H )+^ 

B = Pf - p H 
dans laquefle Pf est la pression de forma- 
tion et p H est la pression hydrostatique de la boue. 



6. Le proc£d£ selon une revendication prgcgdente 
quelconque, dans Iequel la vitesse de variation 25 
dpi de la pression d'admission Pj est d6termin§e 
^ rinstant^et la valeur deC, est prise 6gale k la- 
dite vitesse de variation dp h 



7. Le proc6d6 selon Tune des revendications pr6c6- 30 
dentes quelconque, dans Iequel la vitesse de va- 
riation dp de ia pression de la boue est determi- 
ne k un instant ^ et le volume V D de I'inf lux est 
calcuie par liquation: 

v = Ph^V^p 35 
0 d m gv g -dp 
dans laquelle p H est la pression hydrosta- 
tique de la boue, est la compressibility de la 
boue, V m est le volume de la boue dans le puits 
de forage, d m est la density de la boue, g est I'ao 40 
alteration due k la pesanteur et v g est la vitesse 
moyenne d'6l6vation de ('influx dans le puits de 
forage. 

8. Le proc6d6 selon Tune quelconque des revendi- 45 
cations 1 & 6, dans Iequel la vitesse de variation 

dp de la pression de la boue est d6termin£e k 
Tinstant tc, le volume V 0 de gain de puits est me- 
sur6 et la vitesse moyenne d'61§vation v g de Tin- 
flux dans le puits de forage est calcutee k partir 50 
de liquation: 

V 0 = P H (Qo-Ql) - PhXmVm 
C 2 (Pf-PH) 

dans laquelle p H est la pression hydrosta- 
tique de la boue, p, est la pression de formation, 55 
X m est la compressibility de la boue, V m est le vo- 
lume de la boue dans le puits de forage. 
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